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ABSTRACT 

We have studied in detail the 0.15-15 GHz radio spectrum of the gamma-ray bi¬ 
nary LS 5039 to look for a possible turnover and absorption mechanisms at low fre¬ 
quencies, and to constrain the physical properties of its emission. We have analysed 
two archival VLA monitorings, all the available archival GMRT data and a coor¬ 
dinated quasi-simultaneous observational campaign conducted in 2013 with GMRT 
and WSRT. The data show that the radio emission of LS 5039 is persistent on day, 
week and year timescales, with a variability < 25 % at all frequencies, and no sig¬ 
nature of orbital modulation. The obtained spectra reveal a power-law shape with a 
curvature below 5 GHz and a turnover at ~ 0.5 GHz, which can be reproduced by 
a one-zone model with synchrotron self-absorption plus Razin effect. We obtain a 
coherent picture for a size of the emitting region of ~ 0.85 mas, setting a magnetic 
field of B ~ 20 mG, an electron density of «e ~ 4 x 10^ cm^^ and a mass-loss rate 
of M ~ 5 X 10“® MQyr^'. These values imply a significant mixing of the stellar 
wind with the relativistic plasma outflow from the compact companion. At particu¬ 
lar epochs the Razin effect is negligible, implying changes in the injection and the 
electron density or magnetic field. The Razin effect is reported for first time in a 
gamma-ray binary, giving further support to the young non-accreting pulsar scenario. 

Key words: binaries : close - gamma rays: stars - radio continuum: stars - radiation 
mechanisms: non-thermal - stars: individual: LS 5039 - instrumentation: interferom¬ 
eters. 


1 INTRODUCTION 

Gamma-ray binaries are extreme high-energy systems consist¬ 
ing of a young massive star and a compact object, which ex¬ 
hibit non-thermal emission from radio to very high energy 
gamma rays (> 100 GeV). Their Spectral Energy Distribu¬ 
tion (SED) is dominated by the MeV-GeV photons (see |Dubus| 
|2013| for a recent review). Only a handful number of gamma- 
ray binaries have been discovered up to now: PSR B1259-63 


( [Aharonian et al.|2005bk LS 5039 ( [Aharonian et al.|2005a| ). 


LS I -h 61 303 dAlbert et al. 

|2006k HESS J0632+057 (|Hin- 

ton et al.|2009 

Skilton et al. 

200911, and IFGL J1018.6-5856 

([Fermi LAT Co 

laboration et al.|2012k The first one is the only 


system which hosts a confirmed pulsar. The nature of the com- 
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pact object still remains unknown for the rest of the gamma- 
ray binaries. 

The behaviour of this kind of systems still represents a 
challenge for the high-energy astrophysical community. Al¬ 
though a microquasar scenario, typically seen in High Mass 
X-ray Binaries, was firstly proposed for these sources, the 
most recent results seem to favour a young non-accreting pul¬ 
sar scenario for the four gamma-ray binaries deeply explored 
up to now jDubus|2013[ >. Under the frame of this scenario, the 
broadband non-thermal emission of these extreme energetic 
systems is thought to arise from the particles which are accel¬ 
erated in the shock between the wind of the compact object 
and the wind of the companion star. The gamma-ray emission 
would be produced by Inverse Compton upscattering of stel¬ 
lar photons by the most energetic electrons accelerated in the 
shock between the winds. The same population of electrons 
would produce the synchrotron emission observed at radio 
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frequencies (|Dubus||2006| |Durant et al.|2011t |Bosch-Ramon| 
|20^[20T3) . 

LS 5039 (with lAU name J1826-1449 and coordinates 
aj 20 oo = 18''26’" 15.06% djzooo = -14°50'54.3", |Mold6n et af] 
|2012a) is a gamma-ray binary system composed by a young 
06.5 V star and a compact object of 1-5 Mq, orbiting it ev¬ 
ery 3.9 d in an eccentric orbit with e 0.35 (jCasares et al.| 
|2005 1. The system is located at 2.9 ± 0.8 kpc from the Sun 
i Moldon et al.|2012^ . [Paredes et al.H2000^ proposed for first 
time that LS 5039 was associated with an EGRET gamma- 
ray source. H.E.S.S. detected the binary system in the 0.1- 
4 TeV energy range jAharonian et al.||2005a| >, and later re¬ 
ported variability at those energies along the orbital period 
( [Aharonian et al.|2006t . FermifLPCT confirmed the association 
with the EGRET source by clearly detecting the GeV coun¬ 
terpart, which showed also an orbital modulation ( |Abdo et^T] 
|2009| >. The system has been observed and detected with differ¬ 
ent X-ray satellites, showing moderate X-ray variability and 
a power-law spectrum with low X-ray photoelectric absorp- 


tion consistent with interstellar reddening 

Motch et al. 

1997 

Ribo et al.|1999[[Reig et al.|2003| [Bosch- 

iamon et al. 

2005 


ray flux, correlated with the TeV emission, has been reported 
by |Takahashi et al!] ( |2009| l. Extended X-ray emission at scales 
up to 1-2' has also been detected ( [Durant et al.|20ll) . X-ray 
pulsations have also been searched for, with null results [Rea| 
|etal.|2011| l. 

At radio frequencies, [Marti, Paredes & Rib6| ( |1998[ ( con¬ 
ducted a multifrequency study of LS 5039 with the Very Large 
Array (VLA) from 1.4 to 15 GHz with four observations span¬ 
ning 3 months (one run per month). The radio emission of 
the source has a non-thermal origin, and can be described 
by a power-law with a spectral index of cr = -0.46 ± 0.01 
(where the radio flux density at a given frequency, Sy, scales 
as S'v oc v“). These authors reported a flux density variability 
below 30 % with respect to the mean value.[Rlb6 et al.[ jl999[ ( 
and [Clark et ^ [200 1[ ( presented the results of a radio mon¬ 
itoring campaign conducted on a daily basis during approxi¬ 
mately one year with the Green Bank Interferometer (GBI) at 
2.2 and 8.3 GHz frequencies. The emission was similar to the 
one reported in[Martf et al.[([1998[(, exhibiting a mean spectral 


index of or = 


-0 5 +"-^ 

O.J 0.3 


lRib6 et al.]l999 '. 


Clark et al. 


120011 


constrained any orbital modulation to be < 4 % at 2.2 GHz. 
Searches of pulsed radio emission have been also conducted 
with null results, probably because of strong free-free absorp¬ 
tion by the dense wind of the companion at the position of the 
compact object [McSwain et al.|2011^. 


The radio emission of the source is resolved at mil- 
liarcsecond scales using Very Long Baseline Interferometry 
(VLBI). [Paredes et al.[ ( [2000[ [2002^ found an asymmetrical 
bipolar extended emission on both sides of a bright domi¬ 
nant core, suggesting that LS 5039 was a microquasar. Later 
on, [Ribo et ^ [2008[ reported persistent emission from the 
dominant core and morphological changes of the extended 
emission on day timescales, which were difficult to reconcile 
with a microquasar model. Finally, [Moldon, Ribo & ParedesJ 
( [2012b[ discovered morphological variability modulated with 
the orbital phase, and pointed out that a scenario with a young 
non-accreting pulsar is more plausible to explain the observed 


behaviour. In this case the relativistic wind from the pulsar 
shocks with the stellar wind of the massive companion and a 
cometary tail of accelerated particles forms behind the pulsar 
(see [Mold6n et al.[2012b[ and references therein for further de¬ 
tails). 

At low radio frequencies (< 1 GHz) only few observa¬ 
tions have been published up to now. [Pandey et al.[ ( [2007[ ( 
performed two observations with the Giant Metrewave Ra¬ 
dio Telescope (GMRT) simultaneously at 235 and 610 MHz, 
and reported a power-law spectrum with a spectral index of 
a X -0.8 at these frequencies. On the other hand, [^dambe [ 
[et al.[j2008) conducted additional GMRT observations and re¬ 
ported a turnover at ~ 1 GHz, with a positive spectral index 
of Of a: -(-0.75 at the same frequencies (further discussed in 
[Bhattacharyya et al.|2012[ (. These results would indicate, sur¬ 
prisingly, that the flux variability of LS 5039 increases signif¬ 
icantly below 1 GHz (see a possible explanation in the frame¬ 
work of the microquasar model in [Bosch-Ramon[2009[ . How¬ 
ever, as discussed in this work, the change in the spectrum is 
the result of using different and uncomplete calibration proce¬ 
dures. 

The radio spectrum of gamma-ray binaries is dominated 
by the synchrotron emission, which is expected to be self- 
absorbed at low frequencies. The spectrum could also reveal 
the presence of other absorption mechanisms such as free- 
free absorption or the Razin-Tsytovitch effect (also known as 
Tsytovitch-Eidman-Razin or simply Razin effect, [Homby &[ 
[Williams[1966[ l. This last mechanism, the suppression of syn¬ 
chrotron radiation in a plasma, is well known in solar bursts 
studies and widely reported in colliding wind binaries. These 
last systems probably share with gamma-ray binaries the ori¬ 
gin of their emission: the collision between the winds of the 
two components of a binary system and the presence of shocks 
[Dougherty et al.[2003[ . 

Due to the surprising behaviour of the radio emission 
of LS 5039 at low frequencies obtained with a limited num¬ 
ber of observations, we decided to make an in-depth study of 
the source. We have reduced and analysed two monitorings of 
LS 5039 at high frequencies that we conducted several years 
ago, most (if not all) of the publicly available low-frequency 
radio data, and conducted new coordinated observations to get 
a complete picture of the source behaviour at low and high 
radio frequencies. This has allowed us to monitor the vari¬ 
ability of LS 5039 on different timescales (orbital and long¬ 
term), and determine its spectrum using simultaneous and non- 
simultaneous data. This paper is organised as follows. In the 
next section we present the radio observations analysed in this 
work. In Sect. we detail the data reduction and analysis of 
these data, and the results are presented in Sect. We de¬ 
tail the emission models used to describe the observed spectra 
in Sect. and discuss the observed behaviour in Sect. We 
summarise the obtained results and state the conclusions of 
this work in Sect. [7] 


2 RADIO OBSERVATIONS 

The data presented along this work include two high- 
frequency VLA monitorings, all the available archival low- 
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1998 2002 2006 2010 2014 



Figure 1. Summary of all the data presented in this work in a fre¬ 
quency versus time diagram, including the Modified Julian Date 
(MJD, bottom axis) and the calendar year (top axis). Squares repre¬ 
sent the VLA monitorings in 1998 (blue) and 2002 (red), green cir¬ 
cles correspond to the archival GMRT observations, black markers 
represent our coordinated GMRT (diamonds) and WSRT (pentagons) 
quasi-simultaneously campaign (the open pentagon represents a data 
set that could not be properly calibrated, see text). The radio data pre¬ 
sented here span 15 yr. 

frequency radio observations performed with GMRT, and a 
new coordinated campaign that we conducted in 2013 us¬ 
ing GMRT and the Westerbork Synthesis Radio Telescope 
(WSRT)QT]ln Fig.0 we summarise these observations in a 
frequency versus time diagram. The details of all these obser¬ 
vations, including the used facilities and configurations, the 
observation dates and frequencies, can be found in the first six 
columns of Table[T] 

2.1 VLA monitorings 

Two multi-frequency monitorings have been conducted with 
the VLA: one in 1998 jMarti et al.|1998t and another one in 
2002 (unpublished). In addition to the 2002 data set, we have 
reduced again the 1998 data set to guarantee that all the VLA 
data have been reduced using the same procedures. 

The 1998 VLA monitoring consists on four observations 
conducted between February 11 and May 12 at 1.4, 4.8, 8.5 


' The VLA and ATCA data archives contain several isolated observa¬ 
tions of LS 5039 above 1 GHz which have not been analysed in this 
work, because sporadic flux density measurements are not useful for 
the study of short term and orbital flux density modulations. These 
heterogeneous data sets could be useful in future spectral studies. 

^ We have also analysed two unpublished VLA observations con¬ 
ducted on 2006 December 16 and 18 at 330 MHz (project code 
AM877). However, we only obtained upper limits, because the noise 
is too high to clearly detect the source according to the GMRT data 
presented below. For this reason we will not explicitly include these 
observations and the obtained results in the rest of this work. 

^ We have analysed a LOFAR observation performed in 2011 at 
150 MHz during its commissioning. However, the noise and the qual¬ 
ity of the image at this stage were not accurate enough for our pur¬ 
poses. 


and 15 GFIz using 2 Intermediate Frequency (IF) band pairs 
(RR and LL circular polarizations) of 50 MHz each (project 
code AP357). The amplitude calibrator used was 3C 286. The 
phase calibrators used were 1834-126 at 1.4 GHz, 1820-254 
at 4.8 and 8.5 GHz, and 1911-201 at 15 GHz. 

The 2002 VLA monitoring consists on 16 observations 
conducted between September 23 and October 21 (covering 
around 7 orbital cycles) at the same frequency bands and us¬ 
ing the same calibrators than the previous monitoring (project 
code AP444). 

2.2 Archival GMRT observations 

We have reduced all the existing archival GMRT data of 
LS 5039. These data set includes ten simultaneous observa¬ 
tions at 235 and 610 MHz taken between 2004 and 2008: two 
published in |Pandey et al.| ( f2007| l, one reported in |Godambe| 
|et al.|p008^ , and seven additional unpublished observations. 
There is an additional observation at 154 MHz performed in 
2008. The GMRT data were obtained with a single IF, with 
dual circular polarization at 154 MHz (RR and LL), and sin¬ 
gle circular polarization at 235 MHz (LL) and 610 MHz (RR). 
The 154 and 235 MHz data were taken with a 8 MHz band¬ 
width and the 610 MHz data with 16 MHz (with 64 channels 
at 235 MHz and 128 channels at the other bands). In all these 
observations, 3C 286 and/or 3C 48 were used as amplitude cal¬ 
ibrators and 1822-096 or 1830-360 as phase calibrators. We 
have reduced again the published data to guarantee that all the 
GMRT data have been reduced using the same procedures. 

2.3 Coordinated GMRT-WSRT campaign in 2013 

To study in detail the radio emission and absorption processes 
in LS 5039, we obtained for the first time a quasi-simultaneous 
spectrum of the source at low and high radio frequencies 
with GMRT (project code 24_001) and WSRT (project code 
R13B012) on 2013 July 18-22. 

We observed LS 5039 with GMRT at 154 MHz on July 
18, 20, and 22 in three 8-hr runs, and simultaneously at 235 
and 610 MHz on July 19 and 21 in two 4-hr runs. We con¬ 
ducted WSRT observations at 1.4, 2.3 and 4.8 GHz on July 
19 and 21 in two 10-hr runs, which took place just after the 
GMRT ones at 235 and 610 MHz. The GMRT IFs, chan¬ 
nels and polarizations are the same as the ones described 
in Sect. 12.21 The total bandwidths are 16 MHz at 154 and 
235 MHz, and 32 MHz at 610 MHz (double than in the ear¬ 
lier observations reported in Sect. |2.2[ l. The WSRT data were 
obtained with 8 IFs of 46 channels each, with dual circular po¬ 
larizations at 2.3 GHz (RR and LL) and dual linear polariza¬ 
tions at 1.4 and 4.8 GHz (XX and YY). The total bandwidth 
is 115 MHz at all frequencies. We used 3C 286 and 3C 48 as 
amplitude calibrators for both observatories and 1822-096 as 
phase calibrator for GMRT (no phase calibrators are needed 
for WSRT due to the phase stability). 
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frequency with the l-cr uncertainty (3-(T upper-limits in case of non-detection). 
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3 DATA REDUCTION AND ANALYSIS 


All the radio data have been calibrated and analy sed us¬ 
ing st andard procedur es mainly within AlP j^ Obij^ (Cotton 
|2008|, Parseltongu^ jKettenis et al.||20()^ and SPAIV^ i In- 
|tema et al.|2009^ have also been used mainly to run scripts that 
call AIPS tasks to reduce the data. The raw visibilities were 
loaded into AIPS and flagged removing telescope off-source 
times, instrumental problems or radio frequency interferences 
(RFI). The amount of data flagged due to RFI depends strongly 
on the frequency. Whereas at frequencies above 1.4 GHz we 
do not observe a relevant amount of RFI, at 150 MHz we lose 
an important fraction of the data in some observations, even 
up to 30 %. 

The VLA data have been reduced using amplitude and 
phase calibration steps. Self-calibration on the target source 
was successful at 1.4, 4.8 and 8.5 GHz, but failed for most of 
the 15 GHz observations due to the faintness of the source. For 
this reason, we preferred not to self-calibrate any VLA data of 
LS 5039 to compare the results at all frequencies with exactly 
the same data reduction process. 

The GMRT data have been reduced as follows. First, we 
performed an a priori amplitude and phase calibration using 
an RFl-free central channel of the band. Afterwards, we per¬ 
formed the bandpass calibration. Finally, we calibrated all the 
channels of the band in amplitude and phase considering the 
bandpass calibration. We imaged the target source and con¬ 
ducted several cycles of self-calibration/imaging to correct for 
phase and amplitude errors. Given the huge field of view of 
the GMRT images (few degrees), we considered not only the 
Mt)-plane but the full uvw-space during the imaging process 
( [Thompson] 1999] ). 

In addition, for the GMRT data we have also performed 
a correction of the system temperature, T^ys, for each antenna 
to subtract the contribution of the Galactic diffuse emission, 
which is relevant at low frequencies. This T^ys correction was 
performed using archival self-correlated observations of the 
target source and the calibrators, but also using recent data 
taken at our request after the campaign in 2013. A detailed 
explanation of the system temperature correction process for 
the GMRT data can be found in Appendix The obtained 
Tsys corrections are directly applied to the flux densities of the 
final target images. We note that these corrections imply an 
additional source of uncertainty. 

The WSRT data have been reduced using amplitude and 
phase calibration steps. The solutions for the amplitude cal¬ 
ibrators were directly extrapolated to the target source. The 
data were imaged and self-calibrated a few times. We had cal¬ 
ibration problems at 1.4 GHz, and we could not extract any 
result from the corresponding data. 

During the imaging process we have used a Briggs ro- 


The NRAO Astronomical Image Processing System, 
http://WWW.aips.nrao.edu 
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bustness parameter of zero in all cases. For the VLA and 
GMRT data the synthesized beam had a major axis around 2 
times the minor axis. For the VLA in A configuration we have 
obtained geometrical mean values for the synthesized beam 
between 2.3" (1.4 GHz) and 0.18" (15 GHz). In D configura¬ 
tion these values increase up to 58" at 1.4 GHz. With GMRT 
we have obtained synthesized beams between 18" and 5" (at 
154 MHz and 610 MHz, respectively). The WSRT observa¬ 
tions exhibit an elongated synthesized beam with a major axis 
around 10 times the minor axis as a consequence of the linear 
configuration of the array and the declination of the source. 
For WSRT the synthesized beam was 110" x 10" in PA = 5° 
at 2.3 GHz, and 48" x 5.7" in PA = -1° at 4.8 GHz. 

The measurement of the flux density value of LS 5039 
in each image has been conducted using the tvstat task of 
AIPS, by considering a small region centered around the tar¬ 
get. We also measured the root-mean-square (rms) flux density 
of a region centered on the source, but clearly excluding it, and 
considered this to be the flux density uncertainty. However, in 
some cases the noise around the source is larger than the ob¬ 
tained rms (due to poor i«;-plane coverage, calibration errors, 
etc.). In such cases, we repeated the flux density measurement 
a few times using different region sizes, and considered the 
dispersion of the obtained values as a measure of the flux den¬ 
sity uncertainty. In any case, these uncertainties were never 
above 50 % of the corresponding rms values. The jmfit task 
of AIPS has also been used to double check the flux densities 
ofLS 5039. 

To guarantee the reliability of the measured flux densities, 
we have monitored the flux densities of the amplitude calibra¬ 
tor for the VLA observations and the flux densities of several 
background sources detected in the field of view of LS 5039 
for the GMRT observations. The constancy of the obtained 
VLA flux densities (below 3 % at 1.4 and 15 GHz and below 
0.4 % at 4.8 and 8.5 GHz) and the lack of trends in the GMRT 
flux densities, allow us to ensure that the reported variability 
of LS 5039 (see below) is caused by the source itself, and not 
by additional effects during the observations or the calibration 
process. 


4 RESULTS 

LS 5039 appears as a point-like source in all the VLA, GMRT 
and WSRT images, and we summarise the results for all these 
observations in Table [T] In this section we focus first on the 
VLA monitoring performed in 2002 to study the variability 
of LS 5039 at high frequencies along consecutive orbital cy¬ 
cles. Later on, we move to the GMRT observations to study 
the variability at low frequencies. We continue by reporting 
the high frequency spectra from the VLA monitoring of 2002. 
We then combine all the data to show the full radio spectrum of 
LS 5039 from 0.15 to 15 GHz. Finally, we report on the quasi- 
simultaneous spectra from 0.15 to 5 GHz from the coordinated 
GMRT-WSRT campaign in 2013. 
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Figure 2. Multi-frequency light-curves of LS 5039 as a function of the 
Modified Julian Date (MJD) for the VLA monitoring in 2002. En'or 
bars represent 1-cr uncertainties. 

4.1 Light-curve at high frequencies (> 1 GHz) 

Figurej^shows the multi-frequency light-curves obtained from 
the VLA monitoring in 2002 as a function of time. The flux of 
LS 5039 is persistent, and shows variability along all the ob¬ 
serving period (sj 28 d) on timescales as short as one day. The 
trend of the variability is roughly similar at all frequencies. 
Apart from the daily variability, there is an evolution of the 
flux density on timescales of the order of 10 days. test 
provides a probability of variability of 8-40cr, depending on 
the frequency. Since the flux density values encompass those 
of the 1998 monitoring (see |Martf et al.|1998| and our values 
in Table[TJ and show similar average values at all frequencies, 
we conclude that the emission of LS 5039 is similar in both 
epochs (1998 and 2002). 

In Fig. 1^ we plot the same data as a function of the or¬ 
bital phase. The reduced number of observations and the pres¬ 
ence of some gaps (~0.35-0.55, ~ 0.60-0.70, ~ 0.85-1.00) 
does not allow us to clearly report on any significant orbital 
modulation on top of the reported daily and weekly variabil¬ 
ity. We note that at orbital phase ~ 0.8 there is a large disper¬ 
sion of flux densities at all frequencies. We have also binned 
the multi-frequency light-curves with different bin sizes (0.15, 
0.20, 0.25) and different central bin positions to increase the 
statistics, but we have not seen any significant deviation with 
respect to the mean values or larger variabilities at specific or¬ 
bital phases (not even at phase 0.8). To further search for de¬ 
viations at phase ~ 0.8, we have re-analysed the full GBI data 
at 2.2 and 8.3 GHz, consisting on 284 observations spanning 
340 days ( |Rib6 et al.|1999[ [Clark et al.|2001^ , obtaining null 
results. 


4.2 Light-curve at low frequencies (< 1 GHz) 

Figure [^top shows the light curves of LS 5039 at 235 and 
610 MHz as a function of the orbital phase from all the anal¬ 
ysed GMRT observations: the archival GMRT observations 
(2004-2008) and the two GMRT observations from the co¬ 
ordinated GMRT-WSRT campaign in 2013. The source ex¬ 
hibits a persistent radio emission at both frequencies during 



0 (/> = 3.90603 d, MJDo = 51942.59) 


Figure 3. Multi-frequency light-curves of LS 5039 as a function of the 
orbital phase for the VLA monitoring in 2002. The campaign covers 
around 7 consecutive orbital cycles. Error bars represent l-cr uncer¬ 
tainties. 
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Figure 4. Top: light-curves of LS 5039 as a function of the orbital 
phase at 235 and 610 MHz of the archival GMRT observations from 
2004 to 2008 and the two GMRT observations from the coordinated 
GMRT-WSRT campaign in 2013. Error bars represent l-cr uncertain¬ 
ties (resulting from the combination of the image noise and the T^ys 
correction uncertainties). Upper limits at 235 MHz are plotted at the 
3-cr level. Bottom: spectral indexes derived from the data above. We 
observe a mean spectral index of a k 0.5 between 235 and 610 MHz. 

this 9-year period. At 610 MHz LS 5039 shows clear variabil¬ 
ity (6it). At 235 MHz the variability is not significant (Icr) for 
the most conservative decision of considering the upper-limits 
as detections with flux densities corresponding to 3 times the 
uncertainties. 

These smaller significances in the variability with respect 
to the high frequency ones are produced mainly by the larger 
uncertainties of the flux densities, which are a combination of 
the image noise and the dominant T^y^ correction uncertainties. 
When we apply these corrections, we increase the typical 1- 
3 % uncertainties to 8 % uncertainties in the flux densities at 
610 MHz, and from 6-16 % to 40 % at 235 MHz. As the T^ys 
corrections are scaling factors which are equally applied to all 
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the observations and only depend on the frequency, they only 
affect the absolute values of the flux densities, but not their 
relative differences. This means that we can determine the dis¬ 
persion of these data prior to applying the T^ys corrections. In 
this case, we obtain that the source is variable at 610 MHz and 
235 MHz with a confidence level of 35cr and 6cr, respectively. 
At 610 MHz variability on day timescales is also detected (see 
MID 53391-53392 and MID 56492-56494 in Table[T}. 

Therefore, we see the same behaviour than at high fre¬ 
quencies: variability on timescales as short as one day, with 
a persistent emission along the years, and without clear sig¬ 
nals of being orbitally modulated. Observing the 610 MHz 
data (Fig.|^top), a small orbital modulation seems to emerge, 
with a broad maximum located aX. (j) ~ 0.8-0.2 and a broad 
minimum al (p ~ 0.3-0.7. To guarantee that this variability 
is not due to calibration issues, we have monitored the flux 
densities of up to 7 compact background sources in the field 
of LS 5039. Their flux densities vary in different ways, but 
none of them shows the same trend of variability as LS 5039. 
Whether the orbital modulation suggested above for LS 5039 
is produced by the poor sampling or by an intrinsic modula¬ 
tion of the source will remain unknown until further detailed 
observations are conducted. 

The spectral indexes from Fig. [^bottom provide a mean 
of a = 0.5±0.8. Given the large uncertainties, we can not claim 
for variability at more than a 3-cr level. However, if we do not 
consider the uncertainties related to the T^ys corrections we do 
obtain a significant variability of the spectral index, which is 
exclusively produced by the two observed upper-limits. 

4.3 High frequency spectra (> 1 GHz) 

We show in Fig.j^the spectrum of LS 5039 for each observa¬ 
tion within the VLA monitoring in 2002. The spectra can be 
roughly fitted with a power law, but they show a slight curva¬ 
ture below ~ 5 GHz in most cases. Given this curvature and 
the fact that at 15 GHz we can have a larger dispersion due 
to changing weather conditions, we have derived a power law 
for the remaining frequencies (4.8 and 8.5 GHz) to estimate 
the spectral index a. The obtained values are in the range from 
-0.35 to-1.1, with an average value of ff = -0.57±0.12, com¬ 
patible with previously reported values. No signatures of or¬ 
bital modulation are detected. However, around orbital phase 
0.8 we observe the most extreme spectra, with the steepest and 
flattest ones. 


4.4 Non-simultaneous spectrum 

Figure |^top shows the obtained spectrum (0.15-15 GHz) 
from LS 5039 with all the data presented in this work (VLA 
monitorings in 1998 and 2002, the archival GMRT observa¬ 
tions in 2004-2008, and the coordinated GMRT-WSRT cam¬ 
paign in 2013). All these data cover around 15 years of ob¬ 
servations. The source is always detected above 0.5 GHz. At 
235 MHz the source is detected most of the times, but as al¬ 
ready mentioned there are some upper limits. On 2004 August 
3 (</i = 0.18) the rms of the image at this frequency is much 
higher than the flux density expected for the source and thus 
these data have not been taken into account in what follows. 


Table 2. Average value of the flux density (Sy) of LS 5039 and its 
relative variability (i5Sy ■ Sy*) at each frequency, considering all the 
data presented in this work. values have been considered at each 
frequency. In all cases we report a variability below 25 %. At 2.3 GHz 
we only have two observations to compute the mean value, so the am¬ 
plitude of the variability is not relevant and can not be compared with 
the rest of the data. 


v/ GHz 

^ohs 

5y/ mJy 

dSySd! % 

0.235 

11 

23 ±5 

20 

0.610 

12 

39 ±6 

16 

1.4 

17 

36±5 

15 

2.3 

2 

28.3 ± 1.6 

6 

4.8 

22 

22 ±3 

12 

8.5 

20 

16.8± 1.7 

10 

15 

18 

11±3 

25 


At 154 MHz the source is not detected in any observation, nei¬ 
ther combining all the existing data. The spectrum is clearly 
curved, with a turnover around ~ 0.5 GHz. 

We observe that the emission of LS 5039 along the 
years remains persistent with a moderate variability above the 
turnover frequency (below the turnover the uncertainties be¬ 
come larger). Given this behaviour, we averaged the flux den¬ 
sity of LS 5039 at each frequency, considering the 235 MHz 
upper-limits as detections with flux densities corresponding to 
3 times the uncertainties. We show the average values in Fig.[^ 
bottom and quote them in Table We observe a standard de¬ 
viation below 25 % at all frequencies. The average value and 
the variability reported at 2.3 GHz (from WSRT data) are not 
as accurate as in the other cases because we only have two 
observations at such frequency. 

It is worth noting also that we obtained remarkable dif¬ 
ferences between the GMRT results presented here and those 
published in Pandeyrt^alji 2007| >, [Godambe et al.H2008) and 
[Bhattacharyya et ^ (2012 i. All these differences are related 
to the Tsys corrections of GMRT data (explained in detail in 
Appendix [A^. We fir st realised that |Godambe et ar|j2008^ and 
[Bhattacharyya et ^P012[ > did not take into account these cor¬ 
rections, and thus all the flux densities from their data were 
under-estimated. In fact, the values we obtain before applying 
the Tsys corrections are compatible with the published ones in 
these two papers. In the case of |Pandey et al.| ( |2()07^ , the origin 
of the discrepancy is different at each frequency. At 610 MHz, 
these authors did not apply the mentioned corrections (assum¬ 
ing that the contribution of the Galactic diffuse emission at 
such frequency is small, while we have found that it is rele¬ 
vant). At 235 MHz these authors obtained the Tsys corrections 
using the Haslam approximation, which overestimates the flux 
density in this case (see Appendi x [A| an d Table [AT| >. These 
are the reasons why |Pandey et al.| ( |2007[ l obtained a negative 
spectral index while we obtain a positive one at these low fre¬ 
quencies. 


4.5 Quasi-simultaneous spectra 

To avoid the problem of the non-simultaneity of the data in 
the previous spectrum, we conducted coordinated GMRT and 
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Figure 5. Spectra of LS 5039 from the VLA monitoring in 2002 ordered with the orbital phase. The dashed line represents the power law derived 
from the 4.8 and 8.5 GHz data. Signals of a curved spectrum below ~ 2 GHz are visible in most cases. At orbital phases ~ 0.8 we observe the most 
extreme spectra (the steepest and the flattest ones). 


WSRT observations in 2013 July 19 and 21. The results of 
these observations are shown in Fig. [7] Although the two ob¬ 
tained spectra are qualitatively similar, again with a turnover 
around ~0.5 GHz, we observe differences between them. The 
spectrum from 2013 July 19 exhibits a stronger emission and 
a pure power law behaviour from 0.6 to 5 GHz, while on July 
21 the spectrum is slightly curved. The data at 154 MHz were 
not simultaneous (taken every other day on 2013 July 18, 20 
and 22) and provided upper limits on the flux density values 
of LS 5039. The combination of these three epochs does not 
improve remarkably the final upper limit, which is similar to 
the lowest one in Fig.|^ 


5 MODELING THE LS 5039 SPECTRUM 

We will now explain the different models that can fit the ob¬ 
served spectra of LS 5039 in the 0.15-15 GHz range and we 
will constrain some of the corresponding physical parame¬ 
ters. As discussed previously, the high frequency spectra above 
1 GHz can be roughly fitted with a power law with a negative 


spectral index of cr x: -0.5. This, together with the compact 
radio emission seen at VLBI scales, has led to establishing the 
synchrotron nature of the radio source jParedes et al.|2000l >. In 
addition, we have also unambiguously revealed the presence 
of a low-frequency turnover, which takes place at ~ 0.5 GHz. 

To explain the observed radio spectrum of LS 5039 we 
have considered a simple model. The fact that the extended 
emission is only responsible for a small fraction of the total 
flux density jMoldon et al.|2012b| l allows us to consider as a 
first approximation a one-zone model. The absence of a clear 
orbital modulation in the radio flux density implies that the 
characteristics of the emitting region must be similar at any or¬ 
bital phase from the point of view of the observer. The simplest 
model which verifies this condition is a spherically symmetric 
emitting region. For simplicity, we have also assumed that the 
emitting region is isotropic and homogeneous. In this region 
we consider the presence of a synchrotron emitting plasma. 
The turnover below ~ 0.5 GHz could be produced either by 
synchrotron self-absorption (SSA), free-free absorption (FFA) 
or Razin effect. 

Considering a particle injection following a power-law 
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Figure 6. Top: Non-simultaneous spectrum of LS 5039 obtained with 
all the data presented in this work (VLA 1998 & 2002, GMRT 2004- 
2008, GMRT-WSRT 2013). Error bars represent 1-cr uncertainties. 
Upper limits are plotted at the 3-cr level (one at 235 MHz lies out¬ 
side the figure limits). Bottom: Average flux density value at each 
frequency. EiTor bai's represent the standard deviations. The upper- 
limit at 235 MHz outside the figure limits has not been considered 
in the averaging. This average spectrum has been fitted with three 
different absorption models: free-free absorption (FFA), synchrotron 
self-absoiption (SSA), and SSA plus Razin effect (SSA+Razin). Mean 
square en'ors (within the symbol sizes, not shown in the figure) have 
been used as uncertainties in the fitting process (see text), while the 
data at 2.3 GHz and the upper-limits at 154 MHz have not been con¬ 
sidered. 


(as expected from synchrotron emission) with index p and 
normalization K for the mentioned emitting region, the parti¬ 
cle density distribution is n{E) &E = K E~^ dE. This injection 
produces a flux density emission of the form: 


Sv OC 


aj{y) 

An Ky 




( 1 ) 


according to |Longair]j2011^ , where Q is the solid angle sub¬ 
tended by the source, the emissivity is defined as 


J(v) oc a{p) K B(P+h/2 y-(P-i)/2 (2) 


where a{p) is a function that only depends on the index p, B 
is the magnetic field, I is the linear size of the emitting region 
(or the radius in the spherically symmetric case), and = 
is the absorption coefficient, where 

xf ^ oc K B("+2)/2 ^-(p+4)/2 ( 3 ) 

is the absorption coefficient from the SSA (|Longair||2011|l, 
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Figure 7. Quasi-simultaneous spectra of LS 5039 obtained with the 
coordinated GMRT-WSRT campaign in 2013. The data < 1 GHz were 
taken with GMRT and the data > 1 GHz with WSRT. The green 
squares represent the data of 2013 July 19 (0 » 0.9) and the red circles 
the data of 2013 July 21 (^ ^^ 0.4). The black arrows show the 3-(T up¬ 
per limits from the 154 MHz GMRT data of 2013 July 18, 20 and 22. 
The green dashed line shows a SSA model fit to the July 19 data and 
the red solid line shows a SSA+Razin model fit to the July 21 data. 


where b{p) is a function of /?, and 


(4) 


is fhe absorption coefficient from the FFA ^Rybicki & Light- 1 
|man|1979l l, where is the electron number density of the stel¬ 
lar wind and the wind temperature. Note that with these 
definitions = ;r(f’/rf)^, where d is the distance to the source, 
and p = 1 - 2a, where a is the spectral index. 

As the relativistic synchrotron emitting plasma is sur¬ 
rounded by a non-relativistic plasma from the stellar wind 
of the massive companion, a suppression of the beaming ef¬ 
fect from the synchrotron radiation is expected (the mentioned 
Razin effect). In this case, the emission decreases exponen¬ 
tially at frequencies below a cutoff frequency, vr = 20ne ^ 
This absorption can be approximated by an exponential reduc¬ 
tion in the sync hrotron emission of Eq . 0 with a multiplying 
factor of i Dougherty et al. 2003| l. 

Different combinations of the mentioned absorption 
mechanisms have been considered to explain the ob¬ 
served spectra of LS 5039, namely: SSA, FFA, SSA-l-FFA, 
SSA-l-Razin, FFA-l-Razin, and SSA-l-FFA-l-Razin. We have fit¬ 
ted the data and obtained statistics from the residuals, from 
which we have computed the reduced value using the number 
of degrees of freedom (d.o.f.):;^^ =;^'^/d.o.f. 

In the case of the non-simultaneous average spectrum 
we have taken into account all the data except the values at 
2.3 GHz (given that we only have two observations at such fre¬ 
quency, the inferred mean value is not representative enough) 
and except the upper-limits at 154 MHz. In this case we ob¬ 
serve that any of the considered models fits the data with 
x] « 1, implying that their uncertainties are overestimated. 
The reason is that these uncertainties reflect the variability of 
the source, and thus the dispersion of the data, but not the 
uncertainty of the mean value at each frequency. Therefore, 
we need to consider the mean square errors, dividing the stan¬ 
dard deviation at a given frequency by the square root of the 
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Table 3. Parameters fitted to the averaged data (A.S.) with a 
SSA+Razin model. The data on 2013 July 19 (J19) were fitted with 
a SSA model and the data on 2013 July 21 (J21) with a SSA+Razin 
model. The parameters are: p = \ — 2a, Pi — Pi = 

;b(p+ 2 )/ 2 ^ and vr = 20ne 


Fit 

P 

O 

1 

Pi 

[lO^cmG'P+WZj 

VR 

[lO^Hz] 

A.S. 

2.16 + 0.04 

500 + 800 

3 + 5 

4.1+ 0.2 

J19 

1.867 + 0.014 

3.9 + 0.3 

(2.1 +0.9) X 10® 

- 

J21 

2.24 + 0.08 

200 + 600 

0.4+ 1.7 

4.1+ 0.7 


number of measurements at such frequency. When using the 
mean square errors, we observe that there are three models 
which can explain the spectrum: FFA, SSA and SSA+Razin 
(see Fig. [^bottom, where mean square errors, within the size 
of the symbols, are not plotted). The;yJ from these fits are 2.5, 
4.6 and 1.9, respectively. The SSA+Razin model is thus the 
most accurate model to explain the LS 5039 average emission, 
and predicts a spectral index of ctssA+Razin = -0.58 ± 0.02 for 
the optically thin part of the spectrum. Considering the differ¬ 
ent upper-limits at 154 MHz as possible flux density values at 
that frequency (i.e., a flux density of 3cr ± cr), we observe that 
the SSA+Razin is still the best model to reproduce the data 
(although if the assumed flux density at 154 MHz is as low as 
< 7 mjy, a FFA+Razin seems to be needed). 

In the case of the quasi-simultaneous spectra (Fig.|^, we 
have followed the most conservative choice of considering the 
largest upper-limit at 154 MHz as the possible flux density of 
LS 5039 at such frequency (12±4 mJy). We have fitted the two 
spectra using all the possible models mentioned above, and 
assessed their validity with a test. The SSA+Razin model 
can explain both spectra. However, on 2013 July 19 the con¬ 
tribution of the Razin effect is negligible, and thus we observe 
a pure SSA. The FFA model, in contrast, only explains the 
spectrum on 2013 July 21. The SSA+FFA can explain the two 
quasi-simultaneous spectra, although with a small contribution 
of FFA on July 19, but can not explain the average spectrum 
discussed in the previous paragraph. 

In summary, a SSA+Razin model emerges as the most 
reasonable to explain the spectra of LS 5039, with a different 
contribution of the Razin effect, which might not be present 
at all, depending on the epoch. We note that we only have 1- 
3 degrees of freedom in the fitting process, implying that in 
some cases it is subtle to distinguish between models produc¬ 
ing similar spectra. 

For the discussed models, the underlying physical param¬ 
eters are coupled. Therefore, the number of free parameters 
is lower than the number of the physical ones. For the SSA 
model, we only have three parameters to fit: p. Pi = 
and Pi = K f For the SSA+Razin effect model, we 

have the three previous parameters and an additional one: 
Vr = 20«efi *• We can not decouple these physical parame¬ 
ters without making assumptions on some of them. Therefore, 
we prefer to work directly with the coupled parameters for the 
different spectra, which are quoted in Table 


6 DISCUSSION 

First we will discuss on the physical implications of the mod¬ 
eling of the obtained radio spectra for LS 5039. Second we 
will discuss on the radio variability of the source, making a 
comparison with its multi-wavelength emission and the one 
observed in other gamma-ray binaries. 


6.1 Physical properties from the spectral modeling 

In this work we have reported multi-frequency data of LS 5039 
that has allowed us to obtain a non-simultaneous average spec¬ 
trum in the 0.15-15 GHz range, and two quasi-simultaneous 
spectra in the 0.15-5 GHz range, showing in all the cases 
a turnover at ~ 0.5 GHz. The high-frequency emission is 
roughly explained by a power-law with a spectral index of 
a X -0.5 {p ~ 2), which is typical of non-thermal synchrotron 
radiation, although some indications of curvature below 5 GHz 
are observed. All the observed spectra of LS 5039 can be ex¬ 
plained with the models which are summarised in Table 
A SSA+Razin model explains the average spectrum and the 
2013 July 21 spectrum, showing both similar properties (on 
the parameters p ~ 2.2 and Vr » 0.41 GHz), whereas on 2013 
July 19, the spectrum is explained by a pure SSA model, with 
a negligible contribution of the Razin effect (vr 0). In this 
case we observe a significant lower index p - 1.867 ± 0.014, 
a 1-a deviation with respect the average spectrum. The Pi and 
Pi parameters can only be constrained for the 2013 July 19 
spectrum (see Table]^. Therefore, only in the case of p and Vr 
we can claim for changes between the spectra obtained at dif¬ 
ferent epochs. It seems that changes in the injection index but 
also in the electron density or the magnetic field are present in 
LS 5039 on timescales of about days. 

The Razin effect is widely observed in colliding wind 
binaries (CWB), where there is a non-relativistic shock be¬ 
tween the winds of two massive stars, with typical values of 
Vr ~ 2 GHz jVan Loo|2005t . Interestingly, the vr value in¬ 
ferred for LS 5039 on the average spectrum and on 2013 July 
21 is only a factor 5 smaller. The presence of the Razin ef¬ 
fect in this gamma-ray binary, and the fact that shocks are 
present in CWBs, provides further support to the scenario of 
the young non-accreting pulsar, where shocks also take place, 
for LS 5039. We remark that the presence of the Razin ef¬ 
fect is also supported by the curvature reported on most of the 
LS 5039 spectra below 5 GHz, which is easily explained with 
the consideration of this effect but hardly explained by a SSA 
or FFA model. 

The magnetic field B can be estimated from the Pi pa¬ 
rameter obtained on the 2013 July 19 spectrum, by assum¬ 
ing a solid angle fl for the emitting region, with a radius i in 
the spherically symmetric case considered here. From|Mold6n| 
|et al.|p012b^ it is known that most of the radio emission arises 
from the compact core, which is not resolved with VLBI ob¬ 
servations at 5 GHz, and thus has an angular size with a ra¬ 
dius < 1 mas. At the ~ 3 kpc distance of the source this rep¬ 
resents ~ 3 AU, or about 10 times the semimajor axis of the 
binary system. We consider a solid upper-limit for the angular 
radius of ~ 1.5 mas, and a lower-limit of ~ 0.5 mas 

taking into account the orbital parameters and the absence of 
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Figure 8. Electron number density for the non-relativistic plasma («e) 
as a function of the angular radius of the emitting region {I d~^, bottom 
axis) and of the magnetic field value {B, top axis). B is determined 
from the angular radius and the Pi parameter of the pure SSA fit to 
the 2013 July 19 data (see Table[^. The red dashed line denotes the 
values infen'ed from the vr parameter of the SSA+Razin fits (which is 
directly related to B). Solid black lines represent the values obtained 
with different mass-loss rates M. Values below the red dashed line 
imply an unrealistic mixing above 100 % of the non-relativistic wind 
inside the synchrotron radio emitting relativistic plasma. The vertical 
dotted line shows the central value of 1 mas. 


orbital modulation (which would be produced by the changing 
shock geometry and by absorption effects). From the central 
value of 1 mas we estimate B x 35 mG. For angular radius in 
the range of 0.5-1.5 mas we obtain a wide range of magnetic 
field values of 2-180 mG (see x-axes in Fig.[^. These values 
range from much smaller to close to the B ~ 200 mG value 
inferred from the VLSI images at 5 GFlz assuming equiparti- 
tion between the relativistic electrons and the magnetic field 
( [Paredes e t al.| |2000^. They are also encompassed with the 
3-30 mG range quoted in |Bosch-Ramon| ( |2009| ( using lumi¬ 
nosity arguments for the lower limit and imposing FFA to de¬ 
rive the upper bound. The obtained magnetic field value of 
B X 35 mG for 1 mas allows us to constrain the Lorentz fac¬ 
tor of the electrons emitted at 0.5 GFlz to y ~ 60, confirming 
that we are in the relativistic regime jPacholczyk|1970[ ). Using 
the Vr value obtained from the SSA-l-Razin fits we derive the 
electron density for the non-relativistic plasma, ne. as a func¬ 
tion of B (red dashed line in Fig. [^, with values in the range 
of 4 X 10“*-4 X 10^ cm“^ for the 0.5-1.5 mas range, and of 
7 X 10^ cm“^ for the central value of 1 mas. We can compare 
these values with the ones obtained from the stellar wind ve¬ 
locity and mass-loss rate. The velocity of the stellar wind is 
Li„ X 2440 km s“* (McSwain et al. 2004 1 . Considering a m ass- 
loss rate of M x 5 x 10“^ Mq yr^* Casares et al.|2005 i we 
estimate x 2.6 x 10® cm“^ for 1 mas, and thus a mixing 
of ~ 25 % of the non-relativistic wind inside the synchrotron 
radio emitting relativistic plasma. 

All these values allow u s to estimate the free-free o pac- 
ity, Tff oc M- I Rybicki & Lightman|l979 '. As¬ 

suming a temperature for the stellar wind of x: 1.3 x 10“* K 
(jKrticka & Kubat||2001^ we obtain Tg x 50 for an angular 


Figure 9. Mass-loss rate (M) for different free-free opacities (tr) as 
a function of the angular radius of the emitting region bottom 

axis) and of the magnetic field value (B, top axis). The relation be¬ 
tween and B is the same as in Fig. M Only the region below 
Tff ~ 1 should be considered according to the observations, implying 
a low value for M. The vertical dotted line shows the central value of 
1 mas. 


radius of 1 mas (see Fig. [^, which is not compatible with 
the presence of radio emission. Flowever, from optical spec¬ 
troscopy there are indications of clumping, which might re¬ 
duce the mass-loss rate by one order of magnitude down to 
M X 5 X 10“* Mq yr“* (Casares et al. in preparation). This 
would yield tr 0.5 (Fig.|^, and thus the region would be op¬ 
tically thin to free-free opacity. This lower value of M would 
imply We ~ 2.6 x 10® cm“®, which is smaller than the value 
derived from the SSA-(-Razin model fits (see Fig.[^. For this 
mass-loss rate, only angular radii of 0.85 mas or smaller 
would be supported (or the mixing would be above 100 %). At 
the same time, from Fig.Mwe observe that for sizes < 0.8 mas 
we obtain Tr > 1. Therefore, the inferred value for the radius 
of the emitting region is ~ 0.85 mas, implying a significant 
mixing of the non-relativistic wind inside the synchrotron ra¬ 
dio emitting relativistic plasma, even close to ~ 100 %. An¬ 
other explanation in agreement with this value is the produc¬ 
tion of the radio emission by secondary electrons produced 
after photon-photon absorption within the stellar wind ^Bosch- j 
[Ramon, Khangulyan & Aharonian|2008| l. 

In summary, from the spectral modeling assuming an ho¬ 
mogeneous one-zone emitting region we derive an angular ra¬ 
dius of ~ 0.85 mas (£ ~ 2.5 AU at 3 kpc), yielding B ~ 20 mG, 
and ~ 4 X 10® cm“®. The mass-loss rate must be around 
M ~ 5 X 10“* Mq yr“*, i.e. about one order of magnitude 
smaller than the value reported in [Casares et al.|(2005^ , thus 
supporting the presence of a clumpy stellar wind. A high mix¬ 
ing with the relativistic plasma is supported. 


6.2 Variability 

In this work we have reported day to day variability, trends on 
week timescales and the absence of orbital variability. The ab¬ 
sence of orbital modulation in the radio emission of LS 5039 
puts some constraints on the characteristics of the emitting re- 
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gion or radio core, in the sense that they can not change dra¬ 
matically along the orbit or that the changes in different phys¬ 
ical parameters should compensate each other, always as seen 
from the point of view of the observer. In addition, the reported 
stability along the years of the average radio emission at all 
frequencies implies that this emitting region must be also sta¬ 
ble on these timescales. On the other hand, hydrodynamical 
instabilities in the shocked material or in the outflow prop¬ 
erties could explain more easily the observed variability on 
timescales of days and the trends observed on timescales of 
weeks jBosch-Ramon|200^ . 

In contrast to the radio emission, LS 5039 shows orbitally 
modulated flux at other wavelengths. The X-ray light-curve is 
periodic, exhibiting a maximum around orbital phases 0.7-0.8 
(during the inferior conjunction of the compact object). It also 
exhibits small spikes at 0 0.40 and 0.48 and a strong spike 

?l\. (f) X 0.70, probably originated by geometrical effects. The 
X-ray emission, including these spikes, presents a long-term 
stability ( [Kishishita et al.|20'0^ . The TeV light-curve is simi¬ 
lar to the X-ray one, exhibiting a maximum around phase 0.7, 
a minimum flux at superior conjunction, and a strong spike but 
in this case at phase ~ 0.8 jAharonian et al.|2005a| ). The two 
light-curves are correlated, with the TeV variability ascribed to 
photon-photon pair production and anisotropic inverse Comp¬ 
ton scattering and the X-ray variability to adiabatic expansion 
( [Takahashi et al.||2()09^ . The GeV light-curve, also periodic, 
exhibits the maximum emission around the periastron and su¬ 
perior conjunction (cfi ~ 0.0), with the minimum emission dur¬ 
ing the apastron and inferior conjunction (0 ~ 0.6). There¬ 
fore, the GeV emission is almost in anti-phase with the X-ray 
and TeV emission, and is thought to reflect the change in the 
angle-dependent cross-section of the inverse Compton scatter¬ 
ing process and in the photon density jAbdo et al.|2009^ . 

The absence of orbital modulation in radio is thus a par¬ 
ticular case in the multiwavelength emission of LS 5039. How¬ 
ever, it is interesting to note the deviations from the general 
behaviour observed around orbital phase 0.8 (as it happens at 
X-ray and TeV energies). Although in most cases they are not 
significant enough, we observe an increasing flux density at 
high frequencies around phase 0.8 (Fig.|^, a higher flux den¬ 
sity emission at 610 MHz starting just before phase 0.8 and 
extending until phase 0.2 (Fig. |^, and a large dispersion in 
the high frequency spectral index close to phase 0.8 (Fig.[^. 
Since the radio emission of low energy electrons might arise at 
scales significantly larger than the orbital system size, the fact 
that we see hints of variability during the inferior conjunction 
of the compact object at orbital phase 0.8 might suggest the 
presence of an additional component due to Doppler boosting 
or to a decrease in the absorption at radio frequencies when 
the cometary tail is pointing closer to the observer (as invoked 
by |Dhawan, Mioduszewski & Rupen|2006| to explain the radio 
outbursts of LS I -1-61 303). 

To compare the behaviour of LS 5039 with the one de¬ 
tected in the other known gamma-ray binaries, we summarise 
some of their properties in Table Given the reduced num¬ 
ber of known sources and the large differences in their phys¬ 
ical properties, a statistical comparison between them is im¬ 
possible. However, it is remarkable that LS 5039 is the only 
case in which the radio emission is not orbitally modulated 


in a periodic or almost periodic way. At high energies (from 
X-rays to TeV) all of them, including LS 5039, show a peri¬ 
odic emission. That makes even more special the behaviour 
of the radio emission of this source with respect to the rest 
of gamma-ray binaries. This behaviour could be related to the 
fact that LS 5039 presents the shortest orbital period (3.9 d) 
and the smallest known eccentricity (0.35), which naturally 
imply a higher absorption of the inner part of the radio outflow. 
IFGL JlOl 8.6-5856 can probably be considered as the physi¬ 
cally most similar system to LS 5039, because both have sim¬ 
ilar massive stars (spectral type 06 V versus 06.5 V) and the 
shortest orbital periods (16.6 d versus 3.9 d). Unfortunately, 
the eccentricity of IFGL JlOl8.6-5856 is still unknown, and 
in case of being high it could easily explain the observed or¬ 
bitally modulated radio emission (e.g. due to changes in the 
absorption processes), contrary to the case of LS 5039. 


7 SUMMARY AND CONCLUSIONS 

We have presented a coherent picture of the 0.15-15 GHz 
spectrum of LS 5039, solving the discrepancies of the low fre¬ 
quency data reported in previous publications. We have un¬ 
ambiguously revealed the presence of a curvature in most of 
the spectra below 5 GHz and a persistent turnover that takes 
place at ~ 0.5 GHz. As described in this work, the average 
spectrum of LS 5039 can be approximated by a simple model 
with one-zone emitting region, which can be considered ho¬ 
mogeneous and spherically symmetric, radiating according to 
a synchrotron self-absorption model and exhibiting evidences 
of Razin effect. The Razin effect, reported for first time in a 
gamma-ray binary, explains the mentioned curvature below 
5 GHz. As this effect is commonly observed in colliding wind 
binaries, it is expected to be present in case that the emission 
from gamma-ray binaries arises also from a shock originated 
by the winds of the companion star and the compact object, 
as it happens in the young non-accreting pulsar scenario. We 
observe a certain stability for the parameters of the fits that 
are well constrained between the average spectrum and par¬ 
ticular epochs. However, at other epochs the source shows 
a slightly different behaviour that implies changes in the in¬ 
jection and in the contribution of the Razin effect (and thus 
changes in the electron density or the magnetic field). The 
presence of free-free absorption can not be discarded in some 
of the spectra, although the SSA+Razin model is the only one 
which can explain all the observed spectra. For angular radii 
in the range of 0.5-1.5 mas we derive magnetic field values 
of B ~ 2-180 mG. A coherent picture within the one-zone 
modeling, considering reasonable values of free-free opacity, 
is obtained for an angular radius of 0.85 mas, B ~ 20 mG, 
We ~ 4 X 10^ cm“^, and M ~ 5x 10“* Mq yr“*. These values 
imply a significant mixing of the stellar wind within the rela¬ 
tivistic plasma of the radio outflow. This is the first time that 
a coherent picture of the physical properties of the assumed 
one-zone emitting region is presented for LS 5039, including 
the magnetic field value of the radio emitting plasma. 

We have also shown that the radio emission of LS 5039 
is persistent with a small variability on day, week and year 
timescales. This variability is present in all the explored fre- 
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Table 4. Properties of the known gamma-ray binaries and their emission. For each binary system we list the spectral type of the companion, the 
distance to the system, d, the orbital period Porb, the eccentricity, e, the periastron separation, a(l - e), and the apastron separation, a(l + e). We 
also quote the general behaviour of each binary system: the flux density range at 5 GHz, the spectral index a around this frequency, and the shape 
of the light-cuiwe at radio and other wavelengths. P means that the light-curve is orbitally modulated, S that the emission is persistent, and U that 
the source remains undetected. P* means that the light-curve is orbitally modulated but variations from cycle to cycle are also observed. 


Source 

Spectral 

Type 

d 

[kpc] 

Porb 

[d] 

e 

a{\ - e) 
[AU] 

a(\ + e) 
[AU] 

SsGHz 

[mJy] 

a 

Multiwavelength behaviour 
radio X-ray GeV TeV 

LS 5039<‘> 

06.5 V 

2.9 

3.9 

0.35 

0.09 

0.19 

15-30 

~ -0.5 

s 

p 

p 

p 

1FGLJ10I8.6-5856® 

06 V 

1.9 

16.6 

- 

- 

- 

1-7 

~ 0.0 

p 

p 

p 

p* 

LS I +61 303® 

BOVe 

2.0 

26.5 

0.72 

0.1 

0.7 

20-200 

-0.5-H0.5 

p* 

psis 

p* 

psk 

HESS J0632+057® 

BOVpe 

~ 1.4 

315 

0.83 

0.40 

4.35 

0.1-0.4 

~ -0.6 

p* 

psk 

u 

p* 

PSRB1259-63® 

09.5 Ve 

2.3 

1237 

0.87 

0.93 

13.4 

1-100 

-1.0-0.0 

P>}J 

p 

p 

p 


(l)|Aharonian et al.| 

2005a|,|Casares et al.H2005K Kishishita et al.H2009K|Abdo et al.H2009),|Casares et al.|^2012^,|Zabalza et al.H2013); 

(2)|Napoli et al.H201 1), 

Fermi LAT Collaboration et al.H2012), 

An et al.|12013)|Bordas et al.|120131; (3)|Frail & Hjellming|119911,|Paredes et al.| 

9997), 

Strickman et al.| 

1998|, 

Harrison et al.U2000),|Gregory| 

20021,|Casares et al.|120051,|Aragona et al.|120091,|Acciari et al.|1201 ll,|Hadasch | 

1 et al.| 

2012),|Ackermann et al. 

120131; (4)|Skilton et al.|(20091,|Aragona, McSwain & De Becker|120101,|Bongiomo et al.|(201 ll,|Casares et al.| 


j2()12),|BordaretarH20T^pIiTrt al.H2014); (5)|JoSistonerarffT9^f,|Wang, Johnstoir&Manch^ter]p004^^lNeguerue]TetarH2^1^, 

|Dernbska, Kijak & LewandowsldH2012^,|Chemyakova et al.H2014^. 


quency range (0.15-15 GHz), and the relative variation in the 
flux densities is similar at all frequencies, exhibiting a stan¬ 
dard deviation of ~ 10-25 %. The absence of orbitally mod¬ 
ulated variability constrains the characteristics of the radio 
emitting region. The observed variability would be produced 
by stochastic instabilities in the particle injection or in the 
shocked material, but not by geometrical effects due to the 
orbital motion. Although the persistence of the flux density 
was known at 1-yr timescales, we have extend this knowledge 
up to scales of ~ 15 yr. At orbital phases ~ 0.8 we have de¬ 
tected signatures of an increasing flux density trend at high 
frequencies, of the starting of an enhanced flux density emis¬ 
sion at 610 MHz, and of a large dispersion in the high fre¬ 
quency spectra. Although these signatures are not significant, 
it is notable that they all happen at the same orbital phase when 
there is enhanced X-ray and TeV emission. Additional moni¬ 
toring campaigns, specially at low frequencies but also at high 
frequencies could clarify if this behaviour is the result of a re¬ 
duced number of observations or if it is an intrinsic effect from 
the source due to Doppler boosting or changes in the absorp¬ 
tion mechanisms possibly connected with the X-ray or TeV 
emission. 


LS 5039 remains undetected at 154 MHz with the cumu¬ 
lative of ~ 17 hr of GMRT observations, although according 
to the modeling its flux density should be close to our upper- 
limits. A detection with more sensitive interferometric obser¬ 
vations would improve our knowledge about the absorbed part 
of the spectrum. We expect that the higher sensitivity of the 
Low Frequency Array (LOFAR) would allow to clearly detect 
the source for first time at such frequency. In the future, other 
facilities such us LWA, MWA, SKA-low, together with LO¬ 
FAR, GMRT and the new low-frequency receivers at the VLA, 
will significantly improve the sensitivity and resolution of low 
frequency radio observations in the 10 MHz-1 GHz range, 
thus allowing for detailed studies of absorption in gamma-ray 
binaries. 
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APPENDIX A: Tsys CORRECTIONS IN THE GMRT 
DATA 


The power measured by a radio antenna contains the radio 
emission from a wide number of sources. It is common to sep¬ 
arate this power in two components: P oc + T^y^, where Ta 
is the antenna temperature, which includes the contribution 
from the target sources, and Tsys is the system temperature. 
Tsys includes the contributions from the Galactic diffuse emis¬ 
sion or sky brightness temperature), the receiver noise, 
feed losses, spillover and atmospheric emission, all of them 
considered as the effective receiver temperature, (|Crane &| 
|Napier|1989||Sirothia|200^ : 

Tsys = Tsky + T^c (Al) 


While at high frequencies the contribution of the Galactic dif¬ 
fuse emission, is negligible compared to it becomes 
dominant at low frequencies, thus producing very high sys¬ 
tem temperatures. This is particularly relevant when observing 
close to the Galactic Plane. 

To recover the emission from the target sources, T^ys must 
be determined and removed from the measured power P. For 
single-dish antennas and most interferometers, T^ys can be ac¬ 
curately determined through injected calibration signals and 
is automatically subtracted for each antenna. In the case of 
GMRT, the observations are, in general, conducted with the 
Automatic Gain Controller (AGC) enabled to avoid possi¬ 
ble limitations such as saturation produced by RFI or excess 
power. However, this rescales the gain of the system in pres¬ 
ence of strong emission (such as the one of the Galactic Plane). 
In these kind of observations, the determination of Tsys is not 
accurate enough and it requires to be accounted for during the 
calibration of the radio data. The amplitudes of the complex 
visibilities are proportional to the ratio of the antenna temper¬ 
ature to the system temperature {TJTsyP. In a typical radio 
observation, the amplitude calibrator allows us to recover the 
real flux densities of the sky and the target source. In the case 
of GMRT we will end up with the following visibilities for the 
calibrator (C) and the target source (S): 


Vccv 


T-AC) 

T^yACy 


W oc 


Ta(S) 

TayafS) 


(A2) 


During the calibration process, the amplitudes of the complex 
visibilities are correctly scaled (or calibrated) to flux densities 
by applying a scaling factor G oc TsyAC). The scaling factor 
is determined with the amplitude calibrator, which presents a 
known constant flux density along the time. Therefore, in this 
process the contribution of T^yAC) is evaluated and properly 


removed. Schematically, after the calibration process we re¬ 
cover the calibrated visibilities, V, of the amplitude calibrator: 

Vc = GVc (A3) 

However, for the target source we apply the same scaling fac¬ 
tor: 


Vs =G Vs (A4) 

which will only provide properly calibrated visibilities for the 
target source, if TsyAS) = TsysfC). For target sources at 
high Galactic latitudes and/or at high frequencies the Galac¬ 
tic diffuse emission is negligible, and thus T^ys is similar for C 
and For target sources close to the Galactic Plane and/or 
at low frequencies T^yAS) » T^yAC) (because the amplitude 
calibrator is never located close to the Galactic Plane). In this 
case, we will need to quantify the ratio rsys(S)/7’sys(C), or T^ys 
correction, to obtain the correct flux density values in the tar¬ 
get source field. Since LS 5039 lies very close to the Galactic 
Plane, an accurate correction of T^ys is mandatory to properly 
estimate the flux densities of the source at each frequency. 


AI The Haslam approximation 

A common method to implement this correction is based on 
estimating the sky temperature values, T^i^y, from the measure¬ 
ments made by the all-sky 408 MHz survey, which was con¬ 
ducted with several single-dish radio telescopes (|Haslam et al.| 

[TM2I 1. 

Assuming that the Galactic diffuse emission follows a 
power-law spectrum, we can determine the sky temperature 
at a frequency v for any position in the sky (a, 5) from the 
emission reported at 408 MHz by 

The spectral index for the Galactic diffuse emission is usually 
assumed to be 7 = -2.55 jRoger et al.|1999[ >. The sky bright¬ 
ness at the positions of the amplitude calibrator and the target 
source, convolved with the synthesized beam of our interfer¬ 
ometer, allows us to estimate the ratio of the sky temperatures 
between both positions. From eq. IH} we have seen that T^ys 
can be divided in two terms, with Tsky being dominant at low 
frequencies and/or close to the Galactic Plane. The T^c values 
can be obtained from the tabulated data available in the specifi¬ 
cation documents of GMRT and the Tsky values from eq. ( |A5^ . 
With these data we determine the TsysfSJ/TsysfC) ratio, which 
is the Tsys correction. 

However, this method has some problems. First, it as¬ 
sumes a constant spectral index for the Galactic diffuse emis¬ 
sion across all the sky. Second, it assumes that the response 
of our instrument to the Galactic diffuse emission is the same 
as those of the radio telescopes used for the Haslam survey. 
Third, it does not take into account the T^ys dependence on the 
elevation of the Galactic diffuse emission. 


“ Provided that other effects dependent on the antenna position, such 
as elevation, are considered during the cafibration process. 
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A2 Direct measurement of the self-power for each 
antenna 

A more accurate method to implement the T^ys correction 
should only involve measurements conducted with the same 
radio telescope. In the case of an interferometer, one can obtain 
the power measured at any given frequency and elevation for 
each antenna in the array by self-correlating the correspond¬ 
ing data from that antenna. Although the sky brightness is the 
same for all the antennas, the power that comes from is 
different for each antenna, given its internal origin. Thus, con¬ 
sidering the full array as a collection of isolated single-dish an¬ 
tennas, the problem is reduced to determine the power of each 
antenna separately, the so-called self-power, and then average 
all the self-powers to derive the power of the interferometer. 

The self-correlated data of each antenna is proportional 
to the system temperature. Therefore, the ratio between the 
self-powers for each antenna for the target source and the am¬ 
plitude calibrator is equal to the ratio between the system tem¬ 
peratures at these two positions. This ratio is thus a more accu¬ 
rate measurement of 7’sys(S)/7’sys(C) than the Haslam method 
to estimate the contribution from the Galactic diffuse emission. 

As GMRT can nof produce self-correlafed data at the 
same time than standard correlated data, additional observa¬ 
tions have to be conducted. The self-power for each antenna is 
measured at the position of the flux calibrator and the position 
of the source with the same configuration than in the standard 
correlated observations. To avoid possible dependencies with 
the elevation of the sources, both types of observations should 
be conducted when the sources are at similar elevations. 

The self-correlated data shows stronger RFI than the stan¬ 
dard correlated data, and thus a large number of channels must 
be removed through a dedicated flagging process. The ob¬ 
tained amplitudes in different channels show a large disper¬ 
sion. For that reason, it is recommended to take first the aver¬ 
age of the ratios for each antenna and each channel, to min¬ 
imise the dispersion, and then average all the ratios. The use 
of the median, instead of the mean, is also recommended and 
it has been used along our data reduction process. The derived 
value of the ratio is the correction that must be applied to 
the image of the target source field to obtain the correct flux 
density values. 

The main source of error in this method is the large dis¬ 
persion of all the self-powers measured. This can be reduced 
by accumulating more data obtained in different epochs, be¬ 
cause the Galactic diffuse emission is persistent. 

|Sirothia| j2009|l determined the T^ys corrections for 
GMRT with direct measurements of the full sky at 240 MFlz, 
and compared the obtained results with the Flaslam approx¬ 
imation at that frequency. Although the two methods report 
similar values for the whole map, these authors found a high 
scatter in the comparison of the corrections obtained using 
both methods, with an rms in these differences of ~ 56 %. De¬ 
pendencies on the elevation and diurnal/nocturnal time have 
also been found. 

In the case of our target source, we have found the T^ys 
corrections (C) reported in Table | A l] We observe similar cor¬ 
rections at 610 MHz, but different ones at 235 MHz (of about 
60 %), although the Ch and Csp values are roughly compat- 


Table Al. Tsys corrections (C) determined in this work for the position 
of LS 5039 using the amplitude calibrators 3C 286 and 3C 48 with the 
Haslam approximation (Ch) and with the direct measurement of the 
self-powers (Csp). Uncertainties are reported at l-cr level. In this work 
we have used Csp to correct our data. 


v/ GHz 

Ch 

Csp 

154 

~4.0 

1.9 ±0.2 

235 

~ 3.0 

1.9 ±0.3 

610 

~ 1.6 

1.39 ±0.06 


ible at 3-cr level in both cases. At 154 MHz the corrections 
are very different and clearly incompatible. We have seen that 
the Haslam approximation overestimates the T^ys correction 
for the field of LS 5039 at all frequencies. At 235 MHz we 
have obtained an overestimation by an amount in agreement 
with the rms observed by |Sirothia| ( |2009^ at 240 MHz. As 
the direct measurements of the self-powers provide more re¬ 
liable values of the Tsys corrections, and allow us to compute 
the corresponding uncertainties, we have used Csp to correct 
our GMRT data. 
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